Effective Second Order Susceptibility in Photonic Crystals Composed from 

Centro-Symmetric Materials 
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A technique for obtaining efficient bulk second order susceptibility in non centro-symmetric 
Photonic Crystals (PC) made from centro-symmetric materials is discussed. The effect is based 
on the electric quadrupole effect, strong electromagnetic mode deformation and non homogeneous 
contribution to volume polarization from different parts of the PC. The required symmetry breaking 
is introduced on the macroscale of the PC unit cell. The obtained structural Xstr i^ comparable with 
the second order susceptibility of ordinary non-linear materials. Phase matching can be achieved 
by introducing symmetry modulation (Quasi Phase Matching) during fabrication of the PC. 
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Optical materials with high second order susceptibil- 
ity X ^-re frequently required for both fundamental and 
applied research. By utilizing three wave interactions in 
such media, existing coherent radiation sources can be 
extended to almost the entire optical spectrum. Conse- 
quently the importance of these materials can be com- 
pared with the importance of laser itself. 

The number of efficient second order nonlinear materi- 
als is limited. To possess second order susceptibility the 
material has to be non-centro-symmetric. This imme- 
diately eliminates all amorphous materials and crystals 
from 11 of 32 symmetry classes. The choice between 
remaining materials is also limited due to different con- 
straints: value of X i absorption in required spectrum, 
damage threshold etc. 

Different methods were proposed for development of 
enhanced non-linearity. The first is the search for new 
materials with large non-linear response on molecular or 
molecular arrangement |l| scales. The second is construc- 
tion of composite materials [g| . In this case the effective 
enhancement is achieved due to redistribution of the en- 
ergy inside composite structure. Unfortunately it can im- 
prove only x^'^\ ^'^d not x*-^-*. Recently the use of Pho- 
tonic Band Gap (PBG) crystals made from non-linear 
materials was proposed H , . All these methods require 
symmetry breaking on microscopic atomic or molecular 
dimensions. 

In this Letter we describe a technique for obtaining effi- 
cient second order susceptibility in non centro-symmetric 
Photonic Crystals (PC) made from centro-symmetric 
materials. The effect is based on the electric quadrupole 
transition, strong electromagnetic mode deformation and 
different contributions to the volume polarization from 
different parts of the PC. The required symmetry break- 
ing is introduced on the macroscale of the PC unit cell. 

A local second order polarization P*^^) can be obtained 
in centro-symmetric materials due to the quadrupole ef- 
fect 1^. In quadrupole transitions the required symme- 
try breaking is obtained by asymmetry of the electro- 
magnetic field spatial mode, rather than the asymmetry 
of the electron wave function as in the dipole transition 
case. The second order polarization corresponding to a 
quadrupole transition is: 
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where Q is a fourth-order tensor. Generally integration 
over the volume of such point polarization vanishes, due 
to periodicity of E and gradient dependence of Pq . 
Usually, only some weak signals from interfaces can be 
detected g, |0], however the situation in PCs can be 
quite different. 

PCs are artificial two or three dimensional periodic 
dielectric structures [^, (1). They possess unique opti- 
cal properties for electromagnetic radiation with wave- 
lengths comparable with their period, including the ex- 
istence of full PBGs, anomalously strong dispersion |ic| ] 
and high photon localization near defects fll|. High di- 
electric constant modulation is required to obtain strong 
effects, so several technologies for construction of dielec- 
tric/air 2D (l|, |l| and 3D 0, H, [0 PCs have been 
developed. 

Integration of eq. d3) over the volume in dielectric/air 
PC can be different from zero. The reason is unequal 
contributions to the polarization from different parts of 
the media. The polarization of the air regions can be to- 
tally neglected due to low electron density. Constructing 
PC in such a way that in dielectric part the Pq has one 
sign and in the air the opposite, effective "structural" 
volume polarization can be obtained. 

(2) 

Effective second order susceptibility xh induced by a 
quadrupole effect can be quite large [||: 
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where Xn' ~ \Q\ /^ A is the radiation wavelength, x is 
some characteristic second order susceptibility, d is char- 
acteristic interatomic dimension, I is the modulation scale 
of electromagnetic mode (see Fig. nf) and 77 is some nu- 
merical coefficient. I can be as small as the radiation 
wavelength A and d w 0.3n?7i. The origin of 77 > 10 is 
that the dipole transition matrix element between states 
of different parity always has smaller numerical coeffi- 
cient than quadrupole transition based on the matrix el- 
ement of the same parity Q. The experiments on the 
surface non-linear effect in Si, high index material which 



is used for PCs fabrication, show that in this material 

(2) 

Xn is only about three orders of magnitude weaker than 

X^^) of GaAs for A « l/xm- [0|. It means that properly 
designed PCs made from Si, with non-vanishing contri- 

bution of pQ in all the volume, can artificially provide 
second order susceptibility at least of the same order of 
magnitude as low second order coefficient uniform mate- 
rials (e.g. quartz). The result may be even larger due to 
stronger modulation of the electromagnetic field in PC 
than at an ordinary interface and due to the influence of 
photonic band structure on electron transitions H] . 
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FIG. 1. ID defect lattice, with lattice constant L, in Pho- 
tonic Band Gap (PBG) environment. A single defect in PEG 
environments can possess localized modes for frequencies in- 
side the gap. This means that for the same frequencies in 
periodic lattices of defects the mode can be not localized, but 
strongly modulated. I is the characteristic modulation length. 

A specific realization of the desired material requires 
PCs with highly modulated electromagnetic mode. It 
can be achieved directly in some PCs, e.g. high mode 
modulation was predicted for silicon inverse opal struc- 
ture ijl^, or by introducing periodic defects into the PBG 
substrate. The first approach may be realized using lower 
index materials, because the requirement for PBG does 
not exist. Unfortunately there were no intensive numeri- 
cal search for photonic structures with high mode modu- 
lations and at these moment PCs of this type are not well 
explored. In contrast PBG crystals require high refrac- 
tive index materials, but there are many known struc- 
tures and almost any defect (extra or missing material, 
deviation from the periodicity) possesses locaHzed modes 
for frequencies inside the gap, leading to modulated mode 
structure (see Fig. nh. 

Also in the latter approach the required asymmetry 
can be achieved by proper defect's design. Independent 
of the PBG substrate the defects should be asymmetric 
and possess dielectric/air structure. In the case when 
the electromagnetic mode is omnidirectional in dielectric 
part of the defect, the maximum volume contribution of 



polarization eq. (|l|) can be obtained. In our opinion, 
the hollow cavities, partially filled with the substrate's 
material, are the best candidates for proposed method. 

Let us now consider some specific example (see Fig. H) : 
Optical Parametric Oscillations (OPO) in a lattice of 
ID defects (hollow waveguides) in 2D PBG environment. 
Such waveguides array can be fabricated using semicon- 
ductor lithography technology. To obtain a 2D PBG 
crystal, a hexagonal pattern of deep air columns can be 
etched in the high refractive index substrate, e.g. Si. 
During the same lithographic step the waveguides can be 
introduced by etching parallel stripes. Using quasicrys- 
talline arrangement of air columns instead of hexagonal 
pattern, 2D PBG can be achieved in lower index materi- 
als, e.g. glass |M. 



PBG E Air □ Dielectric 



.(2) 



t 



■l'\ 





" /,co.C()-Ac< 


'm^mii^mm^ 


< 

< 

< 

< 




MJ 


// 






mm^^^m 










^mM'^^i^M^^^^ 


►■ 



FIG. 2. Periodic lattice of waveguides (the size of each 
waveguide is I ~ 2nc/ui) in 2D PBG crystal, and the TE 
electromagnetic mode structure inside it. The mode is modu- 
lated due to the strong decay of electromagnetic field in PBG 
crystals. The photonic crystal (PC) is non centro-symmetric 
due to partially dielectric filling of waveguides. Guiding waves 
in such waveguides is an intrinsic property of PBG crystals. 
The figure is a top view of the proposed 2D photonic crysal 
and can be considered as an image of the lithography mask, 
required for its fabrication. 

Inversion symmetry breaking is achieved by partial 
waveguide filling. It can be done by leaving an unetched 
stripe near each etched one(see Fig. |^). The possibil- 
ity to guide electromagnetic modes in such waveguides 
is intrinsic property of PBG crystals. Guiding in PBG 
crystals has no crucial dependence on the internal waveg- 
uide's structure, as in index guided modes, because the 
confinement of the electromagnetic mode exists only due 
to the external PBG environment. 

Three electromagnetic waves (each TE mode): pump, 
signal and idler, propagate along the waveguides direc- 
tions. Let us assume that pump and signal have almost 
degenerate frequencies uj and lu — Ato inside the forbidden 
gap. The idler's frequency Ato is assumed to be outside 
the gap (see Fig. ||). 

Propagation of some specific radiation in a non-linear 
waveguide, can be described under slow varying ampli- 
tude approximation as pO||: 
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where A is an amplitude {E — AF{p)/y'Dexp{iKz — 
iujt)), to and K are the frequency and wavevector of the 
mode, c is the speed of hght, F (p) is the normahzed 
field distribution in the transverse plane, p is transverse 

coordinate and P (p) is the induced non-linear polar- 
ization. For our specific model (see Fig. g) the equation 
for the signal eq. ^) can be written as: 
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using eq. (1) for F and that the field of TE modes 
has only x direction. 

Under the assumption that the idler wavelength is 
larger than the defect's size, Ex can be assumed to be 
constant in the transverse plane. Under the same as- 
sumption the modes of pump and signal can be taken 

(p) 
to be identical with some maximum E^maxOn the dielec- 
tric/air interface inside the waveguide and zeros on its 
boundaries ||ll|, |gl| (see Fig. g). In this case the integral 
in eq. (0) does not depend on the exact electromagnetic 
mode distribution, but on the field value at the boundary 
points of integration. The integral is tacken only over the 
filled part of the waveguide and can be evaluated analyt- 
ically: 
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For an ordinary x'^^ process eq. (||) becomes: 
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To compare eqs. (|^) and (||) one must take into account 
that the energies of pump modes should be the same in 
both cases: 
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The ratio of effective non-linearities in these processes is: 
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It can be rewritten using eq. (^ ) and that Z « A as: 
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where /^overlap — yQL indicates imperfect modes overlap 
in eq. @}. It can not affect the process seriously, since 
L and I can be of the same order { L k, 21 \s realistic) 
due to strong light confinement in PBG crystals, hence 
previously obtained estimations are valid for this specific 
PC. 

For efficient conversion phase matching has to occur, 
Afc = in eq. (P. Otherwise the signal from different 
points along the propagation will interfere destructively. 
Generally in most materials this condition can be influ- 
enced only by change of the propagation direction. 

Phase matching can be achieved artificially by periodic 
modulation of the sign of the non- linear tensor coefficient. 
It can be achieved in ferroelectric crystals only by peri- 
odic poling. The method is called Quasi Phase Matching 
[ p2[ . In non centro-symmetric PCs modulation of the 
non-linear coefficient can be introduced during fabrica- 
tion (see Fig. pi). 
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FIG. 3. Realization of Quasi Phase Matching in PCs for 
second order non-linear processes, 2-n /A. — Ak. The phase 
matching is achieved by period modulation of the asymme- 
try's direction. It can be done by mask modification during 
lithography step. In 3D the task can be completed by lithog- 
raphy multilayer construction. The desired asymmetrically 
filled cavity can be constructed by changing the defect's im- 
age from negative to positive in different layers. 

These results can be extended to the general case of 
periodic lattice defects inside a PBG Crystal. As it was 
stated earlier single defects in PBG environments can 
possess localized modes for frequencies inside the gap. 
This means that for the same frequencies in periodic 
lattices of defects the mode can be not localized, but 
strongly modulated. Assuming that the electric field near 
the defects is E'max and between them is i?min ~ 0, due 
to eq. (nl) the scaling for structural polarization is: 
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where I is the effective localization length of the mode and 
N is the dimension of the defect's lattice. This result is 
valid for the optimal case of field distribution, e.g. as in 
the previously considered example. The dimension of the 
defect lattice can be smaller than the space dimension. 



The periodicity of the entire crystal impHes translational 
symmetry along all unconfined coordinates. Hence these 
coordinates can be omitted (see Fig. |l|). 

The " structural" polarization has to be compared with 
the ordinary one: 
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where L is the defect lattice constant. The ratio between 
-E-max and ii'moan Can be determined taking into account 
that the total energies of the compared modes has to be 
equal, similarly to eq. (0): 
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where a is some numerical coefficient of order 1. Taking 
into account eq. (0) one can get: 
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The additional numerical coefficient /3ovcriap corre- 
sponds to imperfect overlap of different modes [ p3[ . In 
case of three wave interactions its scaling can be deter- 
mined from eqs. ( p^ , pl| , [T^ ) to be: 



/3. 



overlap 



1 



(14) 



in the case of all waves being confined. Practically 
this can be difficult due to the finite forbidden gap 



in PCs 

\Q\E, 



In case only two waves are confined P^ 



(2) 



inca.nEms,J^ ^ and thc scafiug is: 
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where N is the dimension of the defect's lattice. For 
A^ = 1 results (|l^) and (^ are equivalent. The obtained 
results are not limited to mode modulation by defects 
lattice, but the same formulae can be applied in the case 
of "pure" PC strong mode modulation Q. 

PC fibers [g4| and so called "Super Mirror" |2^] waveg- 
uides are other good candidates for implementation of 
second order non-linear processes based on an electric 
quadrupole transition. They are suited for guiding 
light mode in the air, consequently the partial filling of 
these waveguides can provide effective second order non- 
linearity. Recently very interesting non-linear processes 
were observed in air/silica microstructure optical fibers 
p6[ . Three dimensional defects are also highly interest- 
ing due to the analogy between 3D defects and molecules. 
All these cases require more extensive numerical investi- 
gation. 

The measurement of the effects value can be completed 
by light propagation inside a properly designed PC, e.g. 
in the proposed waveguide's array structure, or by light 
scattering from the surface of PC. In this case bulk and 
surface non-linearities can be distinguished by contribu- 
tion of different light polarizations to the signal. In the 
case of surface non-linearity of centro-symmetric mate- 
rials the tangential component of electromagnetic field 



does not contribute to generated harmonics of incident 
light. 

Integrated optics is one of the possible applications for 
structural Xstr rnaterials. Use of standard non-linear ma- 
terials in integrated optics is difficult. Generally integra- 
tion is either impossible due to incompatibility of dif- 
ferent processes or expensive. The most convenient are 
polymers, but several specific disadvantages (e.g. life- 
time) prevent their broad use. The method described 
above opens new possibilities for incorporating non-linear 
elements in optical chips by e.g. obtaining effective non- 
linearity from the chip substrate itself. 

It was shown that a non centro-symmetric PC, made 
from substrates without bulk second order nonlinearity 
(amorphous materials or centrosymmetric crystals) can 
posses effective volume x comparable with the second 
order susceptibility of ordinary nonlinear materials. The 
effect is based on the electric quadrupole transition and 
local enhancement of WE. The symmetry breaking is 
introduced on the macroscale of the unit cell of the PC, 
instead of atomic scale. The estimations were made for 
A — lfj,m, and for smaller wavelengths the non-linearity 
should be stronger. Quasi Phase Matching can be intro- 
duced in the time of fabrication, so effective non-linear 
processes are feasible. 

Discussions with Dr. G. Berkovich and Prof. Y. Sil- 
berberg are gratefully acknowledged. 
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